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flow in the distribution systems is no longer moving in a single direction and 
this resulted the system to become as active distribution networks (ADN). 
Keywords: One of the main problems in ADN is the voltage regulation issue which is to 
maintain the voltage to be within its permissible limits. Several methods of 
voltage control methods are available and focus is given in finding the 
optimal voltage control using artificial intelligence techniques. This paper 
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Distributed Belek ation presents an optimal and coordinated voltage control method while 
Firefly algorithm minimizing losses and voltage deviation of the network. The optimal and 
Optimal settings coordinated voltage control scheme is implemented on an IEEE 13 bus 


distribution network for loss and voltage deviation minimization in the 
networks. Firefly Algorithm (FA) which is a known heuristic optimization 
technique for finding the optimal solution is used in this work. The results 
are compared with another optimization method known as Backtracking 
Search Algorithm (BSA) for identifying the best setting for solving the 
voltage regulation problem. In order to solve the multi-objective optimization 
issue, the MATPOWER load flow simulation is integrated in the MATLAB 
environment with the optimization algorithm. 
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1, INTRODUCTION 

The large-scale implementation of distributed generation (DG) in distributed systems has resulted in a 
conversion from the inactive arrangement, with unidirectional power flow, towards active distributed systems, 
with miulti-directional power flow. With the presence of distributed generations (DGs), the bidirectional 
scheme has changed from the downstream unidirectional scheme to a more effective and active network that is 
called active distribution networks [1, 2]. A lot of challenges occurred due to the change towards active 
distribution system and one of the key issues that has risen with high dispersion of renewable energy sources 1s 
voltage rise [2]. With uncertain variation of distributed generation (DG) and load, it may lead to larger voltage 
fluctuation thus making the voltage quality issues becoming more serious. Due to the unbalanced nature of the 
renewable energy sources such solar irradiance and wind energy, it can remarkably transform the system’s 
voltage profile thus interfering with the established voltage control scheme, for example, the on-load tap 
changer (OLTC) [3, 4]. The problem of voltage rise is more critical when there is no demand in the system, as 
all the generations are exported back to the distribution networks. Since the objective is to improve the 
negative impact of connecting distributed generation (DG) into the system, essentially there is a need to 
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develop an improved management and control plans to compromise the voltage variation problem in order to 
maintain the voltage not to exceed its permissible range [5, 6]. In order to achieve the objectives of delivering 
good quality of power to the consumers, the task of determining the best possible technique to control voltage 
must be achieved. 

Currently, there are various control plan which have been proposed by previous studies in order to 
overcome the effects of integration of distributed generation in distribution system in terms of voltage 
regulation. Conventionally, a distribution system required a connection of all the distributed generation units to 
operate in constant power factor control (PFC) mode [7-10]. In [10], the research work has proposed two local 
controllers for voltage profile regulations at buses where the renewable energy sources such as wind are 
connected. The first technique used fuzzy logic where it is based on reactive power compensation for 
supporting the distributed generation unit voltages. Another technique relies on the sensitivity analysis where it 
has resolved the problem of under voltages due to the miss-coordination between OLTC and DG local 
controlled. Research works in [11] has addressed the voltage regulation issue in finding the ideal tap setting of 
a voltage regulator by adjusting multiple feeders. Optimal voltage regulation method utilizing active and 
reactive power injection by using distributed generation unit was done in [12]. In another research work [13] 
the optimal power flow is applied in automatic power flow management of England distribution system. The 
outputs of distributed generation are determined by real-time optimal power flow under the thermal constraints. 
Metaheuristics optimization methods such as firefly algorithm (FA) and backtracking search algorithm (BSA) 
have been applied in this research work to find the ideal solutions for coordinated voltage control settings of 
three different voltage control methods in a network integrated with distributed generations. Firefly algorithm 
(FA) which was found in [14, 15] originated from the idea of flashing performance of fireflies which resulted 
in the occurrence of bioluminescent communication. The algorithm operates based on global communications 
among the fireflies. In the research work in [16] it was stated that the firefly algorithm (FA) developed by 
Yang in 2018 is one of the most dynamic method to solve optimization problems in power systems. According 
to the work in [17], FA has been able to obtain and find the overall and local optimal simultaneously and has 
obtained the objectives of meeting the energy load demand required by the consumers with the lowest 
operating costs. Previously the work in [18] demonstrated the ability of connecting capacitor in distribution 
network while achieving loss minimization, voltage improvement and increasing available transfer capability 
(ATC). The research aims at determining the placement of capacitor in 41 bus radial distribution networks by 
using Genetic Algorithm and it has shown that (GA) 1s capable to solve the allocation issues. Based on [19], 
FA is applied in the research work for finding the best ideal solution for time multiplier setting and time of 
operation of relay in the distribution system. Other than that, in [20] the human’s identification techniques 
based on dorsal hand vein where the image of human blood 1s analyzed and Firefly Algorithm [FA] is 
implemented at the stage of identifying people. The identification is carried out by using the distance patterns 
between crossing of fireflies and their matching place. On the other hand, in [21], BSA has been utilized to find 
the optimal size and location of DGs in the distribution system, while reducing the system’s losses and 
improve the reliability of the distribution system. BSA is called a dual-population algorithm because it uses 
current and historical population for solving the optimization problems. It is a simple structure where BSA 
does not rely on the early values of the constraints [22]. 

This paper presents a coordinated and optimal voltage control for voltage regulation using firefly 
algorithm. In this research work, Firefly algorithm (FA) was implemented on an IEEE 13-bus distribution 
network to determine the ideal settings of three different types of voltage control techniques of power factor, 
OLTC and curtailment of power generation by minimizing voltage deviation and power losses. In addition, two 
metaheuristics optimization of FA and BSA methods have been applied to distinguish the achievement 
between the two methods in finding the best optimal coordinated voltage control. The optimization codes have 
been tested in MATLAB in order to evaluate the objective function in reducing power losses and voltage 
variation. 


2. COORDINATED VOLTAGE CONTROL ALGORITHM IN DISTRIBUTION NETWORK 
CONNECTED WITH DGs 

As the capacity of DGs connected in the system increased, the generated power of the DG units may 
surpass the energy consumed by the utilities. The voltage will no longer decreases but somehow it is 
increasing the voltage from the substation along the feeder. It will cause high impact on the line voltage 
because of weak network connection of DGs to the distribution systems. Hence, the controller such OLTC 
and traditional controller may exceed the voltage limits. Therefore, in order to keep the voltage to maintain 
within the appropriate limits, the coordinated voltage control method has been introduced in this research 
work for obtaining the best ideal settings in the power systems. In this section, a coordinated voltage control 
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was developed in this research work for find the ideal voltage control in the power system. The steps of 
developing the optimization codes is presented in 2.2. 


2.1. Determination of Ideal Settings of Coordinated Voltage Control Method 

A constrained nonlinear integer optimization issue has been expressed as a multi-objective 
optimization problem for finding the ideal settings of the OLTC, power factor for distributed generation, and 
the value of active power to be curtailed. The main objective in this work is to reduce the voltage deviation 
and power losses at the load buses. The optimization problem’s fitness function is as follows: 


Finin = YPross) + B Vaev) (1) 


where P,,,; 18 the total power loss, F is the fitness function, y is the coefficient factor for the sum of power 
losses, Vie, 18 the voltage deviation at all system buses, and f is the coefficient factor for Vj.,. The total real 
power losses, P;,,, and voltage deviation, V,,,, are defined using the following equations: 


_— Ploss no DG—Ploss with DG 
Progg = toss no DG" Pioss with DG (2) 


Ploss no DG 


where Pi5s5 no pg 1S the total losses obtained before the connection of distributed generations (DGs) in the 
network, meanwhile P),<5 witn pg 18 the total losses obtained after the distributed generation is integrated in 
the network. 


Pioss = pay Piossist = 152,93) ay -seeg (3) 


where n is the number of lines. 


Viref—Vi 
Vaev = rae (4) 


where V;,¢¢ 1s the indication of voltage at bus 7 with the value used is | and V; is the real voltage at bus 1. 

Hence, according to the power flow equations stated above, the voltage deviation and power losses 
must be improved and reduced. Basically, multi-objective function techniques deliver a bunch of ideal 
results. Therefore, the finest optimization solution for this problem is achieved by using the total of the 
coefficient factor technique. This research work has assumed that the coefficient factor for the overall system 
power losses is 0.6 while the voltage deviation is considered to be 0.4. The values are selected as such with 
the reason that the power loss reduction in distribution network will have a higher economic and technical 
impacts, compared to voltage deviation. 

The inequality constraint is the constraint considered in the optimization problem which involved 
bus voltage in which the bus voltage magnitude must be kept within tolerable operation boundaries and are 
defined as follows: 


Vinin Ss [Vj | = Vinax (5) 


where V,,,;,, is the lower bound of bus voltage limits, V,,,, 1s the upper bound of the voltage limits, and |V;| is 
the root-mean-square value of the i” bus voltage. 


2.2. Implementation of Firefly Algorithm for Finding the Ideal Settings of Coordinated Voltage 
Control Methods 

The optimization using Firefly algorithm (FA) technique has been chosen in this research work for 
solving the ideal settings of the coordinated voltage control methods in the power dispersal network 
connected to distributed generations (DGs). The implementation of Firefly algorithm method has been 
developed and simulated in MATLAB in finding the ideal settings of the coordinated voltage control selected 
[23, 24]. At the same time using the Newton-Raphson load flow, the objective function in decreasing the 
voltage deviation and power losses is simulated in the MATLAB environment. The MATLAB software is 
used for solving the optimization coding which is integrated with the load flow simulation using 
MATPOWER. The flowchart of finding the ideal settings of voltage control in the power distribution system 
is shown in Figure 1. 
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Figure 1. Flowchart of the ideal settings of voltage control with DGs using Firefly Algorithm [24] 


The specific steps for implementation of Firefly Algorithm for solving the ideal settings of the 
coordinated voltage control schemes which are (PFC), (OLTC) and power curtailment are explained as 
follows [25]: 

a) Input the system information such as bus, line and generator data. Initialize the parameters and 
population. 

b) Randomly create the initial positions of distributed generations (DGs) within possible solution 
combination. 

c) Insert the Firefly Algorithm considering the following parameters such as number of fireflies = 50; 
number of dimension = 6; maximum iteration = 1000; alpha = 0.25; gamma = delta = 1. These values of 
alpha, gamma and delta are used in equation stated below for update the movement among all the 
fireflies. 


Xx; =X; + Bo * exp(—yr) . (x; — x;) + ax (rand — =) (6) 


d) Calculate the power loss and voltage deviation by running the MATPOWER load flow 
e) Check for bus voltage magnitude constraint using (5). If the boundary is exceeded, repeat step 4. 
f) Calculate the fitness function using (1). 
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g) Update the bestBrigthness and bestFireflies as an objective function for determining the optimal setting 
problems. 

h) Check the light intensity between firefly i and firefly 7. If the limit is exceeded, continue the step by 
updating the movement of fireflies for finding the best firefly among of the fireflies. These best firefly 
position represent the parameters for OLTC, PFC and generation curtailment. The steps will repeat until 
the maximum population size 1s reached. 

1) The process of optimization is repeated starting from step 4 until the stopping criteria is reached to the 
maximum so as to obtain the best solution. 


3. SIMULATION RESULTS 

The simulation was done by utilizing the IEEE 13-bus distribution test system. The location of 
distributed generations (DGs) are set at buses numbers 645 and 684. The locations of DGs are chosen 
because of bus no 645 and bus 684 have recorded the weakest bus compared to other buses. The load flow 
data was simulated on MATLAB software. The one line diagram of the IEEE 13-bus distribution network is 
shown in Figure 2. In this research work of finding the ideal settings, the number of distributed generation 
units connected in the network is 2. The Firefly Algorithm (FA) method is applied and compared with the 
BSA method for finding the ideal settings of the coordinated voltage control settings of the IEEE 13-bus 
distribution network. 
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Figure 2. Two DGs connected in the IEEE 13-bus distribution networks 


To ensure there is no effect on the voltage profile and losses in the networks, no connection of 
distributed generation is done at the base case system. Only one DG unit is connected to the selected bus 
which are at bus 645 or bus 684 at a time. In this paper, there are 3 types of cases that have been tested and 
the FA and BSA optimization method were coded. These optimizations coded is simulated on the MATLAB 
software. The three cases are explained as follows: 

a. Two DGs with total generation of 1MW are placed in the test system 
b. Two DGs with total generation of 2MW are placed in the test system 
c. Two DGs with total generation of 3MW are placed in the test system 

The objective of determining the voltage deviation, system losses and the ideal values of the 
coordinated voltage control settings which include PFC, OLTC and generation curtailment are recorded and 
compared for both the optimization method applied. 

Table 1 shows the overall optimization results where the values obtained is referred to the power 
factor, on-load tap charger and the generation percentages to be curtailed in the distribution networks. Data in 
Table | for the first case study shows that the settings of tap for FA methods are 1.0299 and 0.9945, whereas 
tap settings for BSA methods are 1.0257 and 0.9830 per unit. Next, the power factor settings for FA are 
0.9063 and 0.9284 meanwhile the settings of power factor for BSA are 0.9015 and 0.9146. The sum of the 
generation percentages to be curtailed using FA is 26% and 23% using BSA. Results for the second case 
study showed that the tap settings using FA methods are 1.0146 and 1.0097 while BSA recorded values of 
1.0289 and 1.0135. The power factor settings using FA are 0.9168 and 0.9397 whereas the values recorded of 
using BSA are 0.9132 and 0.9147. The total percentage to be curtailed using FA is 34% while the percentage 
recorded using BSA is 29%. The third case study also displays consistencies in the achieved results. 
The settings of tap using FA are recorded at 0.9987 and 1.0281 while the ideal tap settings using BSA 
methods are recorded at 1.0300 and 1.0300. The power factor settings obtained using FA, which are 0.9220 
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and 0.9078 compared to BSA method which are obtained the values of 0.9174 and 0.9009. For FA method, 


the total generation percentages to be curtailed is recorded at 28% while BSA method is recorded at 24%. 


Table 1. Optimization Results of FA and BSA for 2 DGs of the IEEE 13-bus distribution system 


2MW 3MW 
FA BSA FA BSA FA BSA 

Best Fitness 0.6627 0.6626 0.6239 0.6234 0.4277 0.4275 
Tap DG1 1.0299 1.0257 1.0146 1.0289 0.9987 1.0300 
Tap DG2 0.9945 0.9830 1.0097 1.0135 1.0281 1.0300 
Pf DG1 0.9063 0.9015 0.9168 0.9132 0.9220 0.9174 
Pf DG2 0.9284 0.9146 0.9397 0.9147 0.9078 0.9009 
Pcurt DG1 14.9% 15.3% 17.11% 19.01% 17.98% 19.56% 
Pcurt DG2 11.64% 8.05% 16.94% 10.27% 10.37% 5.02% 


Based on Figure 3, 4 and 5, the graphs show the convergence characteristics for the three cases with 
IMW, 2MW and 3MW of DGs integrated respectively. The outcome from the simulation showed that the 
BSA method gives the least fitness value compared to FA method and also it has faster convergence than FA. 
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Figure 3. Convergence characteristics of FA and BSA with 2 DGs of IMW generation in the IEEE 13-bus 
distribution system 
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Figure 4. Convergence characteristics of FA and BSA with 2 DGs of 2MW generation in the IEEE 13-bus 
distribution system 
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Figure 5. Convergence characteristics of FA and BSA with 2 DGs of 3MW generation in the IEEE 13-bus 
distribution system 


The main objective of this research work is to improve and reduce the power losses and voltage 
deviation in the distribution systems. Based on Table 2, it clearly explained that after the implementation of 
the optimization method of FA and BSA, the complete performance of the network shows decreases in power 
losses and voltage deviation with optimal distributed generation variables settings. Based on Table 2, the 
power losses for the first case study are reduced to 59.04% using FA and 60.37% using BSA methods. The 
results of power losses for the second case study has reduced to 53.06% and 61.85% using FA and BSA, 
respectively. The third case study has shown a reduction in terms of power losses at 42.68% and 60.37% 
using the optimization methods which are FA and BSA. In terms of voltage deviation, for both optimization 
methods, the percentage of voltage deviation for all cases has recorded to be improved approximately at 11%. 
The percentage calculations of power losses and voltage deviation can be calculated by using the following 
equations: 


Ploss basecase—P loss optimization 
Pic = SS eX 100% (7) 
Ploss basecase 
VDpasecase—VP optimization 
= aco tard hd 0 
VDimprovement = Woes X100% (8) 


Table 2. DG overall results on power loss, voltage deviation, and fitness for all the three cases in 
TEEE 13-bus distribution networks 


DG availability Optimization Method Losses (MW) Voltage deviation (V) Fitness 
Base case - 0.1012 0.9998 - 
with IMW DG BSA 00401 08908 016626 
winamwoa fA =— =< te 
With 3MW DG BSA 0.0401 0.8598 0.4975 


4. CONCLUSION 

As a conclusion, this paper proposed a coordinated optimal voltage control method by using two 
optimization techniques which are Firefly algorithm (FA) and Backtracking Search algorithm (BSA) to 
determine three different voltage control settings, namely, power factor control, on-load tap charger and 
generation curtailment. The two optimizations algorithm selected have been implemented on the [EEE 13- 
bus distribution network connected with distributed generations with different values of generation. The 
results obtained from these two algorithms are compared in terms of the performance of fitness function and 
convergence. From the results obtained, it has been found that BSA performs better compared to FA in terms 
of accuracy by giving lower fitness function value and achieving faster convergence rate. In addition, both 
optimization methods are able to reduce the system losses and also improve voltage deviation in the active 
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dispersal network with the ideal settings of PFC, OLTC and generation curtailment. Thus, all the proposed 
optimization methods have met the objectives of this research work. 
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